Nonequilibrium Extension of the Landau-Lifshitz-Gilbert Equation for Magnetic 

Systems 
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Using the invariant operator method for an effective Hamiltonian including the radiation-spin 
interaction, we describe the quantum theory for magnetization dynamics when the spin system 
evolves nonadiabatically and out of equilibrium, dp/dt =^ 0. It is shown that the vector parameter 
of the invariant operator and the magnetization defined with respect to the density operator, both 
satisfying the quantum Liouville equation, still obey the Landau-Lifshitz-Gilbert equation. 
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I. INTRODUCTION 

Magnetization dynamics has attracted much attention 
both in theoretical and experimental physics. In par- 
ticular, the magnetic information requires a deep and 
fundamental understanding of the dynamics of spin sys- 
tem on a short time scale. The Landau-Lifshitz-Gilbert 
(LLG) equation provides a plausible phenomenological 
model for many experimental results lj. Recently, the 
first three authors have derived the LLG equation from 
an effective Hamiltonian including the radiation-spin in- 
teraction 0- It is assumed there that the magnetic sys- 
tem maintains quasi-adiabatic evolution and obeys the 
condition dp/dt = 0. 

However, for a magnetic system whose Hamiltonian 
H(t) evolves nonadiabatically, its statistics deviates far 
from quasi-equilibrium. Thus the density operator does 
neither satisfy the condition dp/dt — nor is given by 
e -pH(t)_ Instead, the density operator does satisfy the 
quantum Liouville equation 



(1) 



The main purpose of this paper is to derive the LLG 
equation from quantum theory for such a nonequilibrium 
magnetic system with the effective Hamiltonian in Ref. 
0. 

To find the nonadiabatic quantum states for this 
nonequilibrium system, we employ the invariant method 
developed for explicitly time-dependent Hamiltonians in 
the Schrodinger picture One advantage of the in- 
variant operator method is that the eigenstates of the 
invariant operator satisfying the quantum Liouville equa- 
tion are exact quantum states of the time-dependent 
Schrodinger equation up to time-dependent phase fac- 
tors. Further we use this invariant operator to define 
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the density operator Q, in terms of which we define 
the magnetization. Finally we show that the magneti- 
zation satisfies the LLG equation even for nonadiabatic 
and nonequilibrium evolution. 



II. LANDAU-LIFSHITZ EQUATION 

We first consider the spin Hamiltonian without the 
radiation-spin interaction given by 



N 



(2) 



where g is the Lande's g-factor, p,B is the Bohr magne- 
ton, TV is the number of spins in the system, and Si is the 
ith spin operator. The effective field H e // includes the 
exchange field, the anisotropy field, and the demagnetiz- 
ing field as well as the external field. As the Hamiltonian 
has the group SU N (2), we may look for an invariant 
operator with the same group SU N (2) of the form 5] 



N 



l (t) =5^&-Bo(t), 



(3) 



where Ro is a vector parameter to be determined by a 
dynamical equation. The invariant operator, satisfying 
the quantum Liouville equation Q , leads to the equation 



N - ( 
»=i V 



dR 

dt 



+ gH B Ro x H e// =0. (4) 



We thus obtain the equation for the vector parameter 
dR 



dt 



-7R0 X H e // 



(5) 



with 7 = guB- 

We note that the eigenstates of the invariant opera- 
tor @ are exact quantum states, up to time-dependent 



2 



phase factors, of the time-dependent Schrodinger equa- 
tion for the Hamiltonian (J2J. Hence the vector parame- 
ter Ro defines the magnetization Ro = Mo of the system 
during the nonadiabatic evolution. When the Hamilto- 
nian W(t) explicitly depends on time, not only the states 
but also the operators change in the invariant method. 
So any physical quantity defined by 



0(t) = Tr{p(t)6(t)} 
has, in general, the time derivative 



(6) 



.9(9, , . . . 



Tr 



at v 



(7) 



where we use Eq. (QJ in the last step. 

In the first case of no radiation-spin interaction, the 
density operator may be given by £| 



Po(t) 



Zo 



-flZb(t) j Z = Tr{t 



h (8) 



where Tq already satisfied Eq. . Here j3 is the inverse 
temperature. Now we may define the magnetization per 
volume in the general case as 



Mo(t) = -Tr{p Q (t)M }, 



(9) 



where V is the volume of the system and M.® is the mag- 
netic moment operator defined by the external field as 



■^o = = 9^b V Si 

Then it follows from Eq. J7J that 
dM 



(10) 



dt 



-7M0 X H e //. 



(11) 



Note that Eq. 1)1 111 is the Landau-Lifshitz equation. 



III. LANDAU-LIFSHITZ-GILBERT EQUATION 

In Ref. |2j the LLG equation is described by an effec- 
tive theory including the radiation-spin interaction. The 
model Hamiltonian 



H(t) = H (t) + m mt (t), 

with the radiation-spin interaction 



(12) 



N 

'Hint = 9I^b Y,Si- («M 2 H e// -MxH e// ), (13) 
i=l 



is an effective theory in the sense that the Hamiltonian 
involves the magnetization to be determined by the the- 
ory itself. As Hint has the same group structure SU N (2) 
as H.q, we still have an invariant operator of the same 
form 



jV 



I(<) = ^5 r R(t). 



(14) 



Then the quantum Liouville equation (JJJ leads to the 
parameter vector equation 



= -g/isRx (l-AoM 2 )H e// +AMxH e// . (15) 



~~dt 



Noting again that the invariant operator (|14fl determines 
exact eigenstates, we may identify R = M for the nona- 
diabatic evolution. Using M • dM./dt — and solving for 
M x dM/dt, we obtain 

dM dM , . 

— = -7M x H e // +aM x — . (16) 

Thus the spin direction of the system even for the nona- 
diabatic evolution obeys the LLG equation. 

In the nonequilibrium case, using the density operator 



p(t) = ± c -0iW Z = Tr{e-^W}, (17) 
we may define the magnetization as 



M(t) = ±Tr{p(t)M(t)}, 



(18) 



where Ai(t) is the magnetic moment operator 
SH 



M 



N 

gp B J2{( 1 - ^uM 2 )Si + XSi x M j(19) 



<5H e // 

The magnetization can be written as 

M = (1 - AaM 2 )Mi + AMi x M, (20) 

where 



Mi(t) = yTr{pM }. 



(21) 



It should be noted that Mo and Mi are weighted with 
Po and p, respectively. Taking the cross product of Eq. 
(l2*0l with Mi and rearranging M x x (Mi x M), we find 
that M is parallel to Mf. 



M 



(l-AaA-f 2 ) + A 2 (Mi-M) 



1 + AM 2 



Mi. (22) 



Hence the last term in Eq. I|20|l drops out, and we get 
M = (1 - AaAf 2 )Mi. (23) 
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The magnetic moment may depend on time through M. 
From Eq. J7J follows the time derivative of M, which is 
given by 



dM 

dt 



hV 



Tr{p{H,M}} 



1 

V 



(24) 



Note that Eq. I|24(l is quite a general result, valid in all 
circumstances, including nonequilibrium evolution. Fur- 
thermore, it has the same form as Eq. (15) of Ref. @, 
where the static limit dp/dt = was assumed. Evaluat- 
ing Eq. H24|) . we finally obtain the LLG equation 



— = - 7 M x H e// 



aM x — -. (25) 
dt ' 



In deriving Eq. (|2*5|l from Eq. lt24"|) we have used M • 
dM/di = 0, which can be checked self-consistently in 
Eq. (|13) . 

A few comments are in order. The Gilbert damping 
term in Eq. i|25|) has the effect of slowing down the pre- 
cession of magnetization and leads to a final constant 
value. Thus the final state is a thermal state with the 
constant magnetization. Then the vector parameter in 
Eq. (|15l) would have the final value 



R(f = oo) = -gfj, B ((l - AaM 2 )H e// + AM x H 



eff 

The invariant operator settles down to the Hamiltonian 
itself, X(oo) = H(oo). Thus the system, after under- 
going a nonequilibrium evolution, reaches a thermaliza- 
tion process towards the final equilibrium with p(oo) = 
e-^(°°)/Z. 



. The time-evolving magnetization enters the Hamilto- 
nian as a parameter and the Hamiltonian thus provides 
an effective theory. When the magnetization proceeds 
rapidly, the system evolves nonadiabatically and out of 
equilibrium. To treat such a nonequilibrium evolution, 
we employed the invariant method to find the equation 
for magnetization. 

The invariant operator, satisfying the quantum Li- 
ouville equation, provides exact quantum states, up to 
time-dependent phase factors, of the time-dependent 
Schrodinger equation. Due to the group structure of 
the effective Hamiltonian ijE} , we were able to find 
the invariant operator, T(t) in Eq. 1)14(1. whose vec- 
tor parameter defines the magnetization. We further 
used the invariant operator to introduce the density 
operator (|17|) . in terms of which the magnetization, 
M(t) = Ti^e-^W^Wj/Z, is defined for nonequilib- 
rium evolution. We showed that the dynamical equa- 
tion for nonequilibrium magnetization satisfies the same 
equation as for the equilibrium case and, therefore, the 
Landau-Lifshitz-Gilbert equation is valid in all cases of 
time development. The nonequilibrium definition of the 
magnetization in this paper has the following physical 
meaning. After the magnetization reaches a final value, 
the invariant operator reduces to the Hamiltonian itself, 
T(oo) = 7i(oo), and the magnetization with respect to 

the density operator, e _/3:r , is nothing but a thermal en- 
semble average of the magnetic moment operator. 
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